The Long Slow Death of the HBT Puzzle 



Scott Pratt 

Department of Physics and Astronomy, Michigan State University 
East Lansing, Michigan 48824, USA 



Abstract 

At the onset of the RHIC era femtoscopic source sizes inferred from two-particle correlations at 
RHIC defied description with hydrodynamic models. This failure, which became known as the 
HBT puzzle, now appears to be solved. The source of the discrepancy appears to be a conspiracy 
of several factors, each of which contributed to making the evolution of RHIC collisions more 
explosive. These included pre-equilibrium flow, using a stiffer equation of state and adding 
viscosity. 



Hydrodynamics represents the foundation of RHIC modeling. The success of ideal hydrody- 
namics in reproducing elliptic flow and spectra inspired the phrase "perfect liquid", as it would 
appear that the quark-gluon plasma is strongly interacting with perhaps the lowest ratio of vis- 
cosity to entropy of any measured substance. However, these models failed to reproduce femto- 
scopic source sizes inferred from two-pion correlations (a.k.a. HBT measurements, named after 
Hanbury-Brown and Twiss who pioneered similar techniques with photons The hydrody- 
namic models El[3]|4 | also fared much worse than purely microscopic simulations |3]|6]|7l[8), 
which did not incorporate soft regions in the equation of state. The failure of hydrodynam- 
ics, along with the success of what seem to be less plausible models inspired the phrase "HBT 
puzzle". Solving the HBT puzzle would require showing that experimentally inferred source 
sizes can be reproduced within the context of a hydrodynamic picture, and without destroying 
agreement with spectra and elliptic flow measurements. Furthermore, the hydrodynamic picture 
should be interfaced with a cascade to model the breakup stage which is well justified at low 
density, when interactions between hadrons become binary. If unusual or unjustified physics is 
added to the model to explain the puzzle, one might simply be exchanging one puzzle for an- 
other. For example, one can better reproduce the data with a compact Gaussian profile for the 
initial energy density [9|, but such a profile is difficult to motivate. 

The solution to the puzzle appears to be the conspiracy of several shortcomings of the original 
hydrodynamic models. Since none of the shortcomings explains more than half of the discrep- 
ancy, they were often overlooked as searches tended to focus on finding a single simple expla- 
nation of the puzzle. The effects described below are of very different origin, but they all push 
the dynamics in the same direction, increasing the explosivity. More explosive collisions lead to 
more compact phase space distributions of outgoing particles. Since HBT analyses effectively 
measure the size of outgoing phase space clouts, and since these dimensions were over-predicted 
by the hydrodynamic models, increasing the explosivity improves the agreement with data. 

Correlations analyses lead to dimensions of the outgoing phase space distributions. If one 
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looks at particles with asymptotic momenta p, their distribution of relative positions is: 



S P (r) = 



J d i r l d i r 2 /(p,ri,?)/(p,r 2 ,0<5 3 (ri - r 2 ) 



(1) 
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This object is referred to as the "source function", though it actually refers to the outgoing phase 
space density /(p, r, f)- It is related to correlation measurements through the Koonin equation, 



where <f> is the well-understood relative wave function. Rather than fully inverting the measured 
correlation function to determine the S p (r), most analyses fit to Gaussian forms for the source 
function, 



and determine the three radius parameters R ont , R s ^ e and Ri ong , which are all functions of p. 
The longitudinal dimension refers to the beam axis, while the outward axis is in the plane of 
p while being perpendicular to the beam axis. For zero rapidity pairs, or if one boosts along 
the beam axis so that p z = 0, the outward direction is parallel to p. The sideward dimension is 
perpendicular to the other two axes. Particles emitted at much different times tend to separate in 
the outward dimension and lead to large values of R out . The exception is when the second particle 
is emitted further outward a time when the first particle is passing by. Despite this qualifier, the 
ratio /? ut/^side is used to gauge the suddenness of the emission. More explosive collisions lead to 
both smaller R ou t/R s ide and smaller values of R\ ong . For semi-boost-invariant emission, particles 
with zero rapidity come only from regions with collective longitudinal velocities that do not 
exceed the longitudinal thermal velocity v^^n, i- e -> Riong(dv z /dz) ~ v t ^ eim , where dv z /dz is the 
longitudinal velocity gradient. For boost-invariant emission, the velocity gradient equals the 
inverse emission time. Thus, ^i ong is sensitive to the average emission time, which is smaller 
for more explosive collisions. It should be emphasized that even though the experimentally 
determined sources were smaller than hydrodynamic predictions, they are still far larger than the 
initial spatial extent of the reaction zone. Longitudinal sizes approach 10 fm, whereas the initial 
nuclei are compressed to a tenth of a fm, and the sideward dimension are also double the initial 
transverse size. Given that the phase space distribution dimensions refer to only those particles 
with a specific momentum, the overall source sizes are significantly larger than the sizes of the 
sub-regions quoted in HBT analyses. This is perhaps the most direct evidence of the dramatic 
expansion of the interacting fireball in RHIC collisions and underscores the folly of ignoring 
longitudinal transverse expansion. 

After nearly a decade of searching for THE solution to the HBT puzzle, it now appears that 
the explanation comes from the conspiracy of multiple features which were missing or misrep- 
resented in the early models. None of these causes represents more than half of the original 
discrepancy by themselves, but after including all the features (or at least the first three below) 
the models reproduce source sizes within the systematic error of the models (often stated as be- 
ing close to 10%). Each of these effects pushes the resulting source sizes in the same direction. 
To illustrate each improvement, Figure [T]presents extracted source sizes for a model with none 
of the features along with results for each feature added incrementally. The model is a viscous 
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hydrodynamic calculation coupled to a hadronic cascade, which models the evolution once tem- 
peratures fall below 170 MeV. By assuming both boost-invariant and azimuthal symmetries the 
model needs only to consider a one-dimensional radial mesh, which makes for rapid calculations. 
After each feature is added the original energy density is rescaled to match the final experimental 
dNch/drj. The speed of the calculations makes it convenient to analyze the impact of the features 
independently, but unfortunately, elliptic flow analyses are precluded by the reduced dimension- 
ality. As a benchmark (open squares in Fig. [TJ a calculation was performed with none of the 
additional features listed below. It fails in very similar ways, both qualitatively and quantita- 
tively, to previous hydrodynamic calculations missing these features. The improvements below 
are then added incrementally. 

1. Adding pre-equilibrium flow. First generation hydrodynamic models initialized the evo- 
lution at a thermalization time, usually chosen between 0.5 and 1 .0 fm/c. However, flow 
can and should develop even before the matter is thermalized. The importance of pre- 
thermalized acceleration for HBT studies has been demonstrated in several studies during 
the last few years lfT0l l6ll9l. and has been investigated in greater detail in regards to other 
observables ifTTl [T2l [131 . Remarkably, this flow is of similar strength for a wide variety of 
pictures of the initial state 1 14 1. The initial impulse of flow results in particle emission that 
is both earlier and more sudden. The closed red squares in Fig. [T]show how this improves 
the calculation relative to the benchmark calculation, by comparing calculations where the 
hydrodynamic expansion was started at 0.1 fm/c vs. 1.0 fm/c. 

2. The equations of state used in many of the early studies was often too soft. In most cases 
a first order phase transition was assumed. For a first order transition the pressure stays 
fixed for a range of energy densities, whose range defines the latent heat of the transition. 
This is inconsistent with lattice calculations which only show a dip in the speed of sound, 
c 2 s - dP/de. Higher pressures lead to more explosive collisions, and the improvement 
(closed green triangles vs. closed red squares) in Fig. [T]is substantial. The calculations 
used simplified forms for the equation of state. For temperatures below 170 MeV, or 
equivalently for energy densities below e%, the equation of state was that of a hadron gas. 
After that point dP/de = c] was fixed for a range of energy densities, L, chosen to be 
1.6 GeV/fm 3 in the default calculation and 800 MeV/fm 3 for the stiffer case. The speed 
of sound was chosen to be zero in the default calculation, i.e., a first-order transition, and 
chosen to be 0. 1 for the more realistic case. Once the energy density exceeded e/, + L, the 
speed of sound was chosen to be c\ = 0.3, consistent with lattice calculations. 

3. Shear viscosity leads to anisotropic stress energy tensors. Since the expansion is stronger 
along the the beam axis, the transverse pressure is increased relative to the longitudinal 
pressure. This also increases the explosivity llT5l[T6l[T7l . For energy densities above e/, + L 
the shear was raised from \nr\js — to 2. For lower energy densities the shear viscosity 
was set equal to that of a hadron gas, and for the intermediate range, £/,<£< e/, + L, the 
viscosity varied linearly with energy density, effectively as a linear interpolation between e/, 
and e/,+L. A Bulk viscosity was also added, see [ 17], but had little effect on the source sizes 
though it reduced radial flow. The change in the source sizes (filled green triangles to filled 
blue circles in Fig. [TJ is again substantial and along with the two previous improvements 
brings the result with the systematic error of the data. 

4. Using a more compact or more Gaussian initial energy density profiles also increases the 
explosivity. Gaussian profiles [9] lead to nearly perfectly linear velocity profiles, as op- 
posed to Wood-Saxon profiles which have higher density gradients near the surface. Hav- 
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Figure 1: (color online) Gaussian radii for three dimensions: S ou t, Rsidc and R{ ong . Data from STAR (red stars) are poorly 
fit by a model with a first-order phase transition, no pre-thermal flow, and no viscosity (open black squares and solid line). 
Data are better reproduced after including all the features (open circles, solid black line). The incremental improvements 
(open colored symbols and dashed lines) are of similar strength: initial flow (red squares), stiffer equation of state (green 
triangles), viscosity (blue circles). Incorporating a more compact initial profile represents the final feature. 
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Figure 2: (color online) Final emission positions and times for particles with transverse momentum of 300 MeV/c along 
the x (outward) axis. Emission is surface dominated and has a modestly positive correlation between x and t. 

ing the profiles steeper at the surface allows the early emission to gain an advantage over 
the later emission from the core, which leads to larger values of R ou t/R\ on g- Although some 
color glass arguments can lead to more compact profiles, there are no good motivations 
for the large effects of [9| which require Gaussian shapes. However, even modestly more 
compact profiles lead to non-negligible changes in the source radii. The effects of employ- 
ing a color-glass inspired [18] profile, vs. a more standard wounded nucleon picture (open 
circles vs. filled blue circles) is apparent in Fig. [T] 
5. A final improvement to the model comes from enacting a more accurate treatment of the 
two-pion wave function for the relative wave applied to the Koonin equation, Eq. Q. 
Most theoretical treatments have ignored the strong and Coulomb interactions between 
the pions as they are compared to experimental results for which the effect of Coulomb 
interactions is largely divided away with what is known as the Bowler-Sinukov procedure 
|[T9ll20l . A truer comparison involves calculating the theoretical correlation functions with 
the Coulomb and strong interaction, then performing the same Bowler-Sinyukov procedure 
on the theoretical correlations. One then compares the resulting correlation functions to 
what one would obtain from Gaussian sources assuming no such complicated interactions, 
consistent with how the data are treated. The more accurate procedure results in source 
sizes that change only at the level of a percent or two IfTTll . but again, the corrections are in 
the same direction as those listed above. 

A second aspect of the HBT puzzle concerns fits with blast-wave models, which are based on a 
picture of thermal emission from a collectively expanding source, parameterized by a breakup 
temperature T, an outer radius R, a breakup time t, a linearly rising transverse collective velocity 
with a maximum, v max , and an emission duration At. These fits revealed outer radii near a dozen 
fm, with emission confined to within a few fm/c of t = 10 fm/c ETl [22). Even though these 
radii are noticably larger than the HBT radii, which reflect only a subset of the overall size, 
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Figure 3: (color online) Model spectra for pions, kaons and pip and spectra (top to bottom, model represented by lines) 
are compared to PHENIX data [29]. Positive/negative species are represented by circles/squares. 

the parameters suggest an unphysically high breakup density. Combining such densities with 
standard hadronic cross sections suggests mean free paths of 1-2 fm. However, breakup should 
not occur until the mean free path is of order of the system size. Figure [2] shows the outward 
coordinate x and the time r at which emission occured for the final model in Fig. [T] and are 
similar to what was seen in J8J. Points are shown only for particles emitted with momentum 
p x = 300 MeV/c. Emission comes mostly from within a few fm of the surface, and there 
exists a modestly positive correlation between x and f, as the average x for the emission points 
moves outward at approximately a tenth the speed of light. A positive correlation prevents those 
particles produced at later times from being strongly separated from those emitted earlier, which 
leads to smaller outward sizes of the outgoing phase space packet. Other physics elements of the 
breakup are also missing from blast wave pictures, such as the differential cooling and collective 
flow of protons and pions once they lose equilibrium l23l . This emphasizes the importance of 
using realistic dynamical models to compare to femtoscopic data and underscores the limits of 
parametric fits. 

The features included above also affect spectra and elliptic flow observables. Unfortunately, 
the symmetries applied in the model used for this study precludes an elliptic flow analysis. As 
shown in Fig. [3] the model does match experimental spectra despite the added explosivity (al- 
though the baryon yields were significantly off). Some of the explosivity might have been bal- 
anced by the bulk viscosity, which had little effect on the HBT radii. A more detailed study 
should be undertaken to see how the various features listed above impact spectra. In addition 
to spectra, elliptic anisotropies and HBT, the entire host of low-/?, observables should be simul- 
taneously addressed to validate a model. This includes more detailed aspects of femtoscopic 
measurements: radii with respect to the reaction plane ll24ll . correlations of other species, es- 
pecially non-identical particles l25l l26l. and non-Gaussian features of the source function |27|. 
Even if all these data are reproduced, it does not fully validate the model. That would require 
an ambitious statistical analysis of the set of model parameters and assumptions, similar to what 
has been applied to cosmic microwave background analyses [28]. Although these goals require 
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significant effort in the coming years, the current analysis has eliminated any puzzle about fem- 
toscopy for the time being, as the experimental radii appear to be satisfactorily described within 
a rather standard theoretical picture of RHIC collision dynamics. 
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